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Abstract 

Classical Coulomb systems at equilibrium, bounded by a plane dielectric wall, 
are studied. A general two-point charge correlation function is considered. Valid 
for any fixed position of one of the points, a new relation is found between the 
algebraic tail of the correlation function along the wall and the dipole moment of 
that function. The relation is tested first in the weak-coupling (Debye-Hiickel) limit, 
and afterwards, for the special case of a plain hard wall, on the exactly solvable two- 
dimensional two-component plasma at coupling F = 2, and on the two-dimensional 
one-component plasma at an arbitrary even integer T. 
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1 Introduction and Summary 

Near a plane wall impenetrable to the particles (hard wall), the charge correlations of a 
classical (i.e. non-quantum) Coulomb system (plasma, electrolyte, . . . ) at equilibrium 
have special features (see review |I|). On one hand, they have only an algebraic decay 
along the wall (while the bulk charge correlations decay faster than any inverse power 
law), and their asymptotic form obeys a simple sum rule. On the other hand, the charge 
correlation function carries a dipole moment (while, in the bulk, there is no such dipole 
just for symmetry reasons), and this dipole moment obeys another simple sum rule. Some 
relation between algebraic tail and dipole moment is expected: it is the asymmetry of the 
screening cloud of a particle sitting near the wall which induces that long-range tail in the 
charge correlation along the wall. In the present paper, we make this relation quantitative. 

The general classical Coulomb systems under consideration consist of s species a = 
1, . . . , s with the corresponding charges q^, plus perhaps a fixed background of density Uq 
and charge density po- Two cases are of particular interest: the one-component plasma 
(OCP) which corresponds to s = 1 (gi = g), po = —qriQ ^ and the two-component 
plasma (TCP) which corresponds to s = 2 (gi = g, g2 = — g), Po = '^o = 0. The presence 
of a solvent is mimicked by embedding the system in a continuous medium of dielectric 
constant e. The walls are made of a material of dielectric constant ew- The particles 
(and the background if any) interact via the Coulomb potential plus perhaps some short- 
range forces. For a //-dimensional system (in what follows, we will restrict ourselves to 
dimensions i/ = 2,3), the Coulomb potential in vacuum at position r, induced by a unit 
charge at the origin, defined as the solution of the Poisson equation, is 

(1.1) 

where tq is a length scale. It should be remembered that, when z/ = 2, the system, 
with logarithmic interactions, is expected to mimic some general properties of three- 
dimensional Coulomb systems, but does not represent "real" charged particles confined 
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to a plane. For two- and more-component plasmas containing pointlike particles, the 
singularity of v(r) at the origin prevents the thermodynamic stability against the collapse 
of positive-negative pairs of charges: in two dimensions for small enough temperatures, 
in three dimensions for any temperature. In such cases, the above-mentioned short-range 
forces (e.g. hard cores) are needed, without effect on the results of this paper. 

We now define some notations. The microscopic densities of charge and of particles of 
species a are given respectively by 

p(r) =po + Xl?""'°('^)' fic{r) ^^5c,aiK'^- n) (1.2) 

a i 

where i indexes the charged particles. The thermal average at the inverse temperature 
13 — l/{kT) will be denoted by (...). At one-particle level, 

p(r) = (p(r)), nM-{na{r)) (1.3) 

At two-particle level, one introduces the two-body densities 



ap(r, r') = (^Sa,aiS/3,ajS{r- ri)5{r^ - rj) \ (1.4a) 

\i^j I 



and the corresponding Ursell functions 

Uapir, r^) = Uaisir, r^) - na{r)n/3{r^) (1.4b) 
The truncated charge-charge correlation (structure function) reads 

S{r, ¥) = {p{r)p{r^)r = mm) - {Kr))W)) (1-5) 

In the case of the OCP, S takes the form 

S{r, r^) = [C/(r, t^) + n{r)5{r - r^)] (1.6) 

For our purpose, it is useful to introduce "conditional" densities. Let na{r\Q, R) be 
the density of a-particles at point r when there is a charge Q fixed at R. Evidently, if 
13 — 1, . . . , s belongs to the set of charged species forming the plasma, it holds 

na{r\qp, r^)np{r^) = Ua/sir, r') (1.7) 



The excess charge density at point r, due to the presence of the charge Q fixed at R, then 
is 

p^^(r|Q, R) = 9a K(r|Q, R) - n„(r)]| + Q5{r-R) (1.8) 

In terms of p*^^, S* is expressible as follows 

5(ry) = 5:g^n^(r^)p'=^(r|g^y) 

= Y.<ipMr)p''''{r'\qp,r) (1.9) 

Here, we consider a semi-infinite Coulomb system which occupies the half-space x > 
filled with a medium of dielectric constant e; we denote by ythe set of {v ~ 1) coordinates 
normal to x. The plane at x = is a hard wall impenetrable to the particles. It may 
be charged by a uniform surface charge density a. The half-space x < is assumed to 
be filled with a material of dielectric constant ew As a consequence, a particle of charge 
q at the point r = (x, y) has an electric image of charge [(e — evi/)/(e + evy)]? at the 
point r* = (—X, y) 0. Due to invariance with respect to translations along the wall and 
rotations around the x direction, 

S{r, r^) = 5(x,x'; \y- = ^(x',x; ly-i/]) (1.10) 

The cases e^y = oo (ideal conductor wall) and ew = (ideal dielectric wall) are special 
and will not be considered here. For finite ew, the charge structure factor S{x, x'; \y— t/|) 
has several general properties: It obeys a condition of electroneutrality 

/•oo r 

J dx' J dy S{x,x';y) = (1.11) 

The Carnie and Chan generalization to nonuniform fiuids of the second-moment Stillinger- 
Lovett condition |^ results for the present geometry in the dipole sum rule [|, |] 

\ dx dx' dyx'S{x,x';y) = -—- -, u = 2,3 (1.12) 

Jo Jo J Zpn[iy — 1) 

The charge-charge correlations decay slowly along the wall ^ . One expects an asymp- 
totic power-law behavior 

5(x,x';y)^^^^, \y\^oo (1.13) 

\y\ 
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where f{x,x'), which as a function of x or x' has a fast decay away from the wall, obeys 
the sum rule [0, PI 



oo roo 



- dx I dx'f{x,x') = - z. = 2,3 (1.14) 

In this work, we establish a general relation between the structure function S and its 
asymptotic characteristics /. Namely, for any value of x > it is proven that 

\ dx' dyx'S{x,x';y) = — 7r(z/ - 1) / dx' f{x,x'), z/ = 2,3 (1.15) 
Jo J ew Jo 

The Ihs of ( [1.15| ) is a dipole moment (in its expression, due to the electroneutrality prop- 
erty ( |1.11| ), x' can be replaced by x' — x). In contrast with the sum rules ( p..l2| ) and ( |1.14| ), 



the relation ( |1.15|) holds for a given x, without integration over it. However, when both 
sides of ( [1.15| ) are integrated over x from to oo, it is seen that the sum rule (|1.14 ) for 



/ is a direct consequence of the dipole sum rule ( |1.12D , and vice versa. 

More generally, it can be assumed that the excess charge density ( |1.8| ) also obeys a 
condition of electroneutrality similar to ( |1 . 1 1|) and has an asymptotic behavior similar to 



(1.13) 

p-(r|g,ig)^ ^^^'^;^\ lyl^oo (1.16) 

\y\ 

where we have chosen R = {X, 0), with F as a function of x or X having a fast decay 
away from the wall. Under these assumptions, we derive the relation 

[ drxp^^'irlQ^R) = —7i{u -1) r dx F{x\Q,X), z/ = 2, 3 (1.17) 
J ew Jo 

valid for an arbitrary Q. Here too, in the Ihs of ( 1.17|) , x can be replaced hj x — X. On 



account of (|1.9| ), this more general relation immediately leads to (|1.15|) with f{x,X) = 
Ef3qfsnf3{X)F{x\qf3,X). 

The paper is organized as follows. Section 2 is devoted to a general derivation of 
the basic result ( |1.17| ). The validity of this result is checked in the Debye-Hiickel limit 



/3 — (section 3) and in two dimensions at the special value of the coupling constant 
r = = 2, for both the TCP (section 4) and the OCP (section 5) in contact 



with a plain hard wall {ew = e = 1). In the case of the two-dimensional OCP, we were 
able to document the validity of formula ( p..l5| ) even for an arbitrary even integer F by a 
non-trivial application of sum rules derived in ref. [0| using the technique of Grassmann 
variables. 



2 General Derivation 

In this section, the abbreviated notation /9^^(r) will be used for the excess charge distri- 
bution p'^^{r\Q, R) defined in ( |1.8|) , with R = {X,0). For simplicity, we first consider 



the case when there are no dielectric media (e = e^y = 1). Our derivation of (|1.171 ) is 



based on the assumption that there are good screening properties in the bulk: the charge 
distribution p'^^(r) is localized near the wall and the electric potential 0(r) it creates in 
the plasma at a macroscopic distance from the wall vanishes. This electric potential is 



J dr^ v{r- r^)p^''{r^) (2.1) 



where v{r — r') is the Coulomb interaction (|1 . 1| ) . When x is large enough for the charge 
distribution p^^{r) to be negligible, v{r — can be expanded in powers of x' (with the 
notation r = (x, y), etc. . . ), and one obtains 

ct>{r) = J dy' v{x, y - y') dx' p^\x\ y') - J dy'^^^^^^^^^ dx' x' p'\x\ y') 
+ ... (2.2) 

where the higher-order terms of the expansion involve higher-order derivatives of v with 
respect to x. 



We shall now consider the Fourier transform of (2^) with respect to y, using the 
convolution theorem. In (|2.2|) , let us define 

a'''{y')= / dx'p^^(x',y') (2.3) 
Jo 

(from a macroscopic point of view, a'^^{'i/) can be regarded as a surface charge density). 
The total charge / dy (j'^^{y) vanishes, as required by the perfect screening of Q. Also, it 
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results from (|1.16| ) that cr^^{y) has the asymptotic behavior 



^^^(y) ^ ^ (2.4) 

where 

PCX) 

A= / dxF{x\Q,X) (2.5) 

JO 

Therefore, the Fourier transform of <j'^^{y) with respect to y has the small wave number 
behavior 

a^^'il) = Jdy exp(-iZ- t/)a"^(2/) ^ -(z/ - l)7rA|Z| (2.6) 

The Fourier transform of v{x, y) is (z/ — l)7rexp(— |Z|x)/|/| and thus the Fourier transform 
of dv/dx is —{v — l)7rexp(— |Z|a;) and so on for higher-order derivatives of v. Using these 
transforms in the convolution theorem gives the Fourier transform 0(x, t) of (|2.2| ), which, 
at Z = 0, is found to be 

0(x, 0) = (z/ - l)7r [-(z/ - l)7rA + P] (2.7) 

where 

P = j dra;p""(r) (2.8) 
Since x has been assumed to be large, 0(a;, 0) = and ( |2.7| ) results into 



(z/ - l)7rA + P = (2.9) 



On account of ( p.5|) and (|2.8| ), ( p.9|) is (|1.17| ) in the special case (e = = !)• 



The generalization to other values of the dielectric constants is straightforward. In 
(p. 21), one must add to p'^^^x',]/) its electric image [(e — ew)/{^ + ^w)]p'^^{~x' , l/)- 
This is equivalent to keeping the integration range (0, oo) for x', but multiplying the first 
integral in (|2.2|) by 2e/ (e + ew) and the second integral by 2evy/(e + ew)- This gives the 
general form of ( |1 . 1 7| ) . 
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3 Debye-Hiickel Limit 

We check the relation (1.17) for the general Coulomb system defined in the Introduction, 
along a plane wall, in the weak coupling limit (3 —>■ (Debye-Hiickel limit). It is convenient 
to introduce the Fourier transform with respect to y of the excess charge density ( |1.8| ) 

p-(x,X, = Jdy exp(-iZ. y)p-(r|g,i2) (3.1) 

One defines a bulk inverse Debye length n by 

f,' = 27i{u-l)(3{J2qln^)/e (3.2) 

a 

where is the bulk density of species a. A minor generalization of the calculation in 
ref. gives pq as a sum of its bulk Debye-Hiickel form plus a "refiected" term: 



e{K' + - ew\l\ 
e(«:2 + /2)i/2 + eH^|/|^^P 



-iK^ + l^y/^ix + X)]^ + Q6ix-X) (3.3) 

The asymptotic behavior of p'^^{y) is determined by the kink of its Fourier transform 
p'^^{l) at |/| = 0, Q{ew/^) Gxp[—hi{x + X)]\l\ (in the sense of distributions, the inverse 
Fourier transform of \l\ is — l/[(t^ — l)7r|y|'^] ). One does find the asymptotic behavior 
( CT ) with 

F{x\Q, X) = - exp[-K{x + X)] (3.4) 

Tc^u — ije 

Therefore 

r dx F{x\Q, X) = , '^^^^ exp(-/tX) (3.5) 

JO 7r(z/ — l)Ke 

For finding the dipole moment associated to p'^^{r\Q, R), one first notes that the 
integral of p*^^ over y is the Fourier transform ( p.3D taken at 1=0. The integral over x is 
easily computed, and one checks that ( |1.17|) is obeyed. 



4 Two-Dimensional Two- Component Plasma 

We now check ( |1.17| ) on the two-dimensional TCP at the special value of the coupling 
constant F = /9g^ = 2, in the case of a plain rectilinear hard wall. Without loss of 
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generality, we shall take the charges ±g as ±1. The corresponding correlation functions 
are known [IC]. We would like to consider the excess charge density (|1.8|) in the special 
case Q = q = 1, i.e. when the particle fixed at i? = {X,0) is one of the particles (say 
a positive one) of the system. Although, at F = 2, for a given fugacity, the densities 
diverge, the Ursell functions are finite. Therefore, instead of p'^^{r\ + 1,R), we consider 
the finite quantity proportional to its non-self part 

n+(i2)[p^^(r| + 1,R)- Q5ir-R)] = U++{r, R) - U-+{r, R) (4.1) 

{Usis2 was called pi^^^ in ref. fl^). In that same reference, the possible surface charge 
density carried by the wall was chosen as —a and we shall keep this choice in the present 
section. 

The model has a rescaled fugacity m [which has the dimension of an inverse length such 
that the bulk correlation length is l/(2m)]. The Ursell functions in (|4.1|) are expressible 
in terms of auxiliary functions gs+{x, X, y), where s = ±, asQ 

U,+ ir,R) = -sm^\gs4x,X,y)\^ (4.2) 

The g functions are given by their Fourier transforms 

/oo 

dy exp{-ily)gs+{x,X,y) (4.3) 
-oo 

which are 

g++{x,X,l) = — I exp (-A;|x - X| ) - exp[-A;(x + X)]|, /<0 (4.4a) 
^++(^,^,0 = g| exp {-k\x -X\) + l~[^fj^ exp[-fc(x + X)]|, / > Qf4.4b) 
g-+{x,X,l) = — |[Z — 27rcr + sign(a; — X)] exp (— /c|s — X| ) 

-(fc + /-27ra)exp[-A;(x + X)]|, /<0 (4.4c) 
g^^{x,X,l) = — 1[/ — 27rcr + A; sign(x — X)] exp (— A;|x — X| ) 

+ (A;-/ + 27ra)exp[-A;(a; + X)]|, />0 (4.4d) 



^We use the symmetry relations (2.15) of ref. the first one of which is misprinted: its rhs should 
be replaced by its complex conjugate. 
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where k = [m^ + {I - 27rcr)^]^/^. 

The asymptotic behavior of the g{y) functions is governed by the discontinuity of 
their Fourier transforms g{l) at / = 0. In the sense of distributions, the inverse Fourier 
transform of sign (Z) is i/^ny). Thus, 

g++{x,X,y) ~ — — — — exp[-ko{x + X)] (4.5a) 

Ztt^Kq — zncrjy 

g^+{x,X,y) ~ -^exp[-A;o(x + X)] (4.5b) 
Iny 



where ko = [m? + (27rcr)2]^/2. Using (|0| ) in and (|]T]) gives 



where 



Therefore 



f/++(r,i^)-f/_+(r,i^)^^^^^4^, lyl^oo (4.6) 

y 



F{x\X) = -p-{ko + 2na) exp[-2fco(x + X)] (4.7) 



/ dx F{x\X) = -—^(ko + 27ia) exp{-2koX) (4.8) 
Jo 47r^ 

For computing the dipole moment 

roo roo 

P= dx{x-X) dy [U++{r, R) - U^+{r, R)] (4.9) 

Jo J-oo 

one first considers the integral over y. With regard to ( [4.2|) and 



/oo jY)'^ roo ^ ^ 

dy [U+^ir, R) - U^^ir, R)] = -— dl {[^++(x, X, /)]' - [^^+(x, X, OjH 
-oo ZTT J-oo 

(4.10) 

In the special case a = 0, using ( |4.4| ) in ( |4.1CI| ) gives for the dipole moment, after some 
calculation, 

jji"^ roo roo ( / ] 

P = / dxix-X) d/ <^exp(-2A;|x-X|) + exp[-2A;(x + X)] , a = 

2tx Jo Jo y m J 

(4.11) 

Performing the integral over x before the one over / gives 

Tfl 

P = exp(-2mX), a = (4.12) 

4:71 
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Comparing ( |4.8| ) when a = with ([4.12|) , one checks that ( |1.17|) is obeyed (here e = ew = 
1)- 

For checking ( [1.17^ in the general case a 7^ 0, it is enough (and easier) to check that 
the derivatives with respect to 27rcr of both its sides are equaL From ( [4 .81) one obtains 

1 / 2na 



d 1 / zttctN 

1(^,1 ^(-1^) = + 1-) (1 -4--^)exp(-2«=oX) (4.13) 

On the other hand, the integrand in the rhs of (|4.1CI| ) depends on / and a only through the 



combination V = I ~ 27rcr, and its analytic form changes at /' = — 27rcr from some function 
/i(Z') to some other function f2{l')i as apparent in ( |4.4| ). Therefore, the corresponding 
integral can be rewritten in the form 

di'Mi')+ di' un (4.14) 

-00 J — 2ttu 

and one obtains 

dJ 



d(27ra) =-/i(2^^) + /2(2vro^) (4-15) 
This gives d/d(2'7rcr) of the integral over y in ( ^.9| ). Then, one computes the integral over 
with the result that dP/d(27r cr) is tt times ( [4.13|) . This completes the check of (|1.17|) 
on the present model. 



5 Two-dimensional one-component plasma 

We consider the two-dimensional OCP in contact with a plain hard wall {ew = e = 1), 
localized at a; = and charged by a "surface" charge density go", first at the coupling 
constant F = fiq^ = 2, then for any even integer F (cr = 0). We work in units such that 
ttuq = 1, where Uq is the background density. We intend to check the relation ( p..l5| ). 



5.1 r = 2 

The one-body density at a distance x from the wall is given by (see eq. (2.16) of ref. 



2 /-oo exp 



(t - xV^f 



n{x) = nQ—= / — — ^dt (5.1 



y/n J-7TaV2 l + $(t) 
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where 



2 1'°° 

^t) = / exp{-u^)du (5.2) 

VTT Jo 

is the error function. The two-body Ursell function (called p^p in eq. (2.18) of ref. 0) is 



U{x,x'; \y - y'\) 



-nl exp 



(X — X 



l\2 



2 roo 6Xp 
^/tI J-ttctV2 



t-{x + x')/V2\ -it{y-y')V2 



dt 



(5.3) 



The structure function S is expressed in terms of n and U in formula ( |1.6|) . The asymptotic 
/-function takes the form (see eq. (2.21) of ref. 0) 



f{x,x') = -n^q 



2 2 2 exp [-2(x + T^af - 2{x' + Trcr)^] 



TT 



It is easy to show that 



— / dx f{x,x') 
Jo 



1 + ^{-^la^/2) 



1 exp [— 2(x + TTcr) 



(5.4) 



^0 ■ ■ y27r3/2 1 + $(-7rav^) 

On the other hand, performing first the integration over the y-coordinate, one obtains 



(5.5) 



— / da;' / dy x'S{x,x';y) = X 



q^ Jo 



it-xV2f 



where 



7r3/2 J-^aV2 1 + $(t) 

8 roo exp 2t^ + 2\p2tx 

exp (-2x ) [l + $(t)]2 



(■OO 

/(t) = / dx' x'exp(-2x'2 + 2y2te') 

JO 



■dt 



I{t)dt (5.6) 



i + [1 + $(t)] exp(t2) 



After simple algebra, one arrives at 

— / dx' \ dy x'S{x, x'\ y) = Ji + J2 
q^ Jo J -00 



where 



Ji 
J2 



I foo I Q 

dtz ^T-T— exp 



72773/2 1 + $(t) dt 



dt 



exp(— 



V^7r2 J-naV2 [1 + $(t)]2 
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exp 



'-{t-x^/2f 
-{t - xV2f 



(5.7) 

(5.8) 

(5.9a) 
(5.9b) 



Using the equality 

exp(— t^) 0r d 
[1 + ~ Ydi [1 + <^{t) 

in J2, the consequent integration per partes imphes 



(5.10) 



J. = - 1 exp[-2(x + ..)^]_ 

y27r3/2 l + $(-7rav^) ^ ^ 



Inserting J2 into (|5l8| ) and comparing with ( ^.5|) one gets the expected relation (|1.15| ) with 
evi/ = e and = 2. 

5.2 r = even integer 

Let the OCP with logarithmic interactions be confined to a compact two-dimensional 
domain V . The positively oriented contour enclosing the domain V , denoted by dV , 
is defined parametrically as x = X{ip),y = Y{ip);(pQ < f < ipi. For example, the 
circle enclosing the disk of radius R centered at the origin admits the parametrization 
X{(p) = Rcosip,Y{ip) = Rsiinp;0 < < 27r. Integrals over the l^-domain can be 
expressed in terms of the (9\^-contour integrals according to the rule 

'dQ dP 



where 



/ P{x, y)dx = r P[X{^),Y{^)]X'{^) 
Q{x,y)dy= rdy,Q[X{y,),Y{^)]Y'i^) 

dV J ipo 



(5.12b) 



The neutralizing background of uniform density uq creates the one-particle potential 
—qnoVo{r) where 

vo(r) = I d^r' v{\r- r'l) (5.13) 
Jv 

The corresponding electric field is —quQEQ^r) where 

Eo{r) = E^,{r)x+El{r)y=-Vv,{r) (5.14) 
with X and y being unit vectors in the x and y directions. For the half-plane of interest, 

VQ{r) = const - vrx^; E^ = 27ix, E^ = (5.15) 
13 



For a disk of radius R centered at the origin 0, 



fo(r) = const — 7rr^/2; = -nx, Eq = iry 



(5.16) 



By mapping the two-dimensional OCP onto a discrete one- dimensional Grassmann 
field theory for the coupling constant F = even integer, two new kinds of sum rules for 
the structure function S were established in ref. [|ll| for an arbitrarily shaped ^-domain. 
The first sum rule reads (see formulae (61a,b) of ref. |ll| where U is called n^) 



-Pn.l,VESir')Sir,r') ^ ^ + //V -f/(.y) 

dy' 



-f3no I dV ^^(r')S(r, r') 
Jv 



(5.17a) 
(5.17b) 



With the aid of the prescription ( |5.12|) , the last terms on the rhs of these sum rules can 
be rewritten as follows 



dV ^U{r, r^) 
V oy 



d<^f/[r; {X,Y)]X' 



(5.18a) 
(5.18b) 



The second sum rule reads (see formula (45a) of ref. ||TT 



-f3n^ [ dV [r' ■Eoir')]Sir,r') 

= 2n(r) + r-Vn(r)+ / U[r; {X,Y)]{XY' - X'Y) (5.19) 

For the half-plane, since the Ursell function U{r, r^) goes to zero at large |r— r^|, eq. 
( FlTap yields 



— 27r/5no / dx' / dy x' S(x,x']y) = — — [ 
Jo J-oo dx J- 



dy U{x,0;y) 



(5.20) 



This is an OCP generalization of the WLMB equations [0, [|T3| which were originally 
derived for neutral systems. It is trivial to obtain the dipole sum rule (|1.12|) by integrating 
both sides of ( ^.201) over x from to cxo, then taking into account that lim^^oo n{x) = hq, 
and finally considering the electroneutrality condition ( |1.11| ) at x = 0. 
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For a disk of radius R, the addition of eq. ( ^.17a|) multiplied by x and eq. ( ^.17b|) 
multiplied by y, with the substitutions (|5.18|) , leads to 

dn(r) 



-n(3n^ I dV (rV)5(ry) =r^^+ r d^' {r-R')Uir,R' 

Jdisk dr JO 



(5.21) 



where, in polar coordinates, R' = {R,ip'). Eq. ( p.l9| ) takes the form 

- n(3no ( dV |rf 5(r, r') = 2n{r) + + i?^ d(^' f/(r, i?' 

Jdisk dr JO 



(5.22) 



Since |r— r^p = |rp + Ir'p — 2r-r', the combination of relations ( ^.21| ) and (|5.22| ) with 
the electroneutrality condition ( |1 . 1 1]) results in 



J disk 



2n(r) — r 



dn{r) 
dr 



r d^' i2r- R! - R^)U{r,R:) (5.23) 

JO 



Let us move the origin to the disk boundary by introducing x = R—r. Eq. ( ^.25 ) becomes 



Tipno r d^' dx'iR-x')\r-r^\'^Siry) 

JO JO 



2n{x) + {R-x 



dn{x) 



2tt 



dif' \2R{R - x) cos{^' -v)- R^] U{r, R'){5.24) 



dx JO 

We now take the limit i? — oo. The Ihs of ( |5.24|) is dominated by the large values of 
|r— r^l. The asymptotic behaviour ( |1.13| ) is equivalent to 



r— r'\ S{r, r') ~ f{x,x') + . . . 



(5.25) 



where the higher-order terms vanish for large |r— /|. Substituting ip' ^ y = R{f' — f) 
on the rhs of (|5.24| ) and collecting all terms of order R, one finally arrives at the equation 

dn{x) 



271 (3nQ / dx' f{x^ x') 



dx 



dy U{x,0;y) 



(5.26) 



formulated for the half-plane. This equation is exactly of the form ( 5.20 ), with the 
expected identification (|1.15|) for the case under consideration e = ew = ^ and u = 2. 
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6 Conclusion 



We have derived a general relation between the asymptotic behavior of a charge correlation 
function along a plane wall and the dipole moment of that correlation function, in the 
general form ( |1 . 1 7| ) . In the particular case of the charge structure factor, that relation 
becomes ( |1.15| ). 

In those exactly solvable models for which the two-body correlations are known, finding 
their asymptotic behavior is often easier than directly computing the dipole moment. 
However, in the present paper, we did compute this dipole moment, for the sake of 
checking the general relation. 



The relation ( 1.15 ) was first observed in the special case of a two-dimensional one- 
component plasma, as one more application of a mapping onto a Grassmann field theory 



TT|. Afterwards, we realized that the relation could be derived in general for any Coulomb 
system. 
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